Hexagons, Kinks, and Disorder in Oscillated Granular Layers 
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Experiments on vertically oscillated granular layers in an evacuated container reveal a sequence of well-defined pattern 
bifurcations as the container acceleration is increased. Period doublings of the layer center of mass motion and a parametric 
wave instability interact to produce hexagons and more complicated patterns composed of distinct spatial domains of different 
relative phase separated by kinks (phase discontinuities). Above a critical acceleration, the layer becomes disordered in both 
space and time. 



PACS numbers: 83.70.Fn, 47.54.+r, 83.10.Pp, 83.10.Ji 



The transport, mixing, and segregation of granular 
materials is important in industries ranging from food 
to mineral processing, yet a basic understanding of the 
physical mechanisms underlying the collective dynamics 
of grains is lacking. Recent experiments on vertically 
vibrated granular materials show a variety of phenom- 
ena including heap formation and convection size 
segregation ||, and traveling waves 0. These phenom- 
ena, although of practical importance, are caused by sur- 
rounding gas and/or sidewall driving Q, leaving open the 
question of whether there is any intrinsic self-organized 
behavior in these systems. We report here robust pat- 
terns that arise spontaneously, not from interstitial gas 
or sidewall forcing, but from correlations induced by mul- 
tiple collisions between the grains and by the coherent 
motion of the particle layer and the container. 

Our experiments on vertically vibrated granular lay- 
ers yield spatial patterns composed of standing waves 
that oscillate at either one-half or one-quarter of the 
drive frequency /. Spatial domains of different relative 
phase separated by phase discontinuities (kinks) appear 
in all patterns except those just beyond the initial in- 
stability. Examples of stripes, squares, hexagons, kinks, 
and a disordered state are presented in Fig. [l| Figure 
U shows the stability region for each pattern as a func- 
tion of / and the dimcnsionlcss acceleration amplitude 
r = 4-7T 2 f 2 A I g , where 2A is the peak-to-peak amplitude 
of the sinusoidal displacement of the container and g is 
the gravitational acceleration. The transitions are well- 
defined, only weakly dependent on /, and non-hysteretic 
(except for the transitions to squares). We will now de- 
scribe the experimental methods and then show how the 
patterns and their associated thresholds arise from the in- 
teraction of two distinct mechanisms: a parametric wave 
instability and period doublings. 




FIG. 1. Patterns in a 1.2 mm deep layer at / = 67 Hz: (a) 
f/2 stripes (r = 3.3), (b) //2 hexagons (r = 4.0), (c) flat 
with kink (F = 5.8), (d) competing f/4 squares and stripes 
(r = 6.0), (e) f/4 hexagons (r = 7.4), and (f) disorder (r = 
8.5). 
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FIG. 2. Stability diagram showing transitions in a 1.2 mm 
deep layer. The vertical dashed lines indicate the frequencies 
above which only stripes appear in the wave regimes. Closed 
(open) square and circular symbols denote the transitions to 
waves with increasing (decreasing) V. 

A layer of 0.15-0.18 mm diameter bronze spheres is 
placed in the bottom of a cylindrical container that has 
inner diameter 127 mm and height 90 mm; the wall and 
lid are Lucite while the base is aluminum to reduce elec- 
trostatic effects. The layer is 7 particles deep except for 
the data in Fig. ||, where a thicker layer (12 particles 
deep) improves the signal-to-noise ratio for the acceler- 
ation measurements. The container is evacuated to 0.1 
Torr, a value at which the volumetric effects of the gas 
are negligible || and heaping is not observed. An elec- 
tromechanical vibration exciter drives the container and 
the resulting acceleration is measured to a resolution of 
O.Olg. Patterns are illuminated at low angle with a strobe 
light and are recorded with a video camera located above 
the container. 

The instabilities leading to the different patterns can 
be understood in terms of two dimensionless parameters 
that characterize the dynamics of the layer: r = ft&t, the 
layer free-flight time, and 7 = v c / gt^ t , which is approx- 
imately the acceleration of the granular layer relative to 
the plate during the time of collision (v c is the relative 
collision velocity) B. We directly measure r, while we 
deduce 7 from a simple one-dimensional (ID) model of 
an inelastic ball on an oscillating plate. Figures ||(a) and 
^(b) show the dependence of r and 7 on the driving ac- 
celeration r, and Table || indicates how the patterns can 
be characterized by the values of r and 7/7^5 where 
7 cr it marks the onset of parametric waves. 

In our model system, the ball is completely inelastic 
because the collisions of the particle layer and container 
are almost completely inelastic due to multiple internal 
grain collisions. The ball motion is computed by as- 
suming that free-flights begin whenever the ball and plate 
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FIG. 3. Comparison of observations for a particle layer to 
calculations for a model consisting of a ID completely inelas- 
tic ball, (a) Flight times for the layer (1.9 mm deep, / = 67 
Hz), T cxp t (A), and for a completely inelastic ball, r ca ic (°). 
(b) The dimensionless layer collision acceleration y/y C rit (the 
average value is displayed where 7 is multivalued), (c) Con- 
tainer acceleration (collisions removed), (d) and (e) Impact 
acceleration for stripes (r = 3.3) and hexagons (r = 4.2) 
respectively (vertical arrows indicate successive collisions of 
regions with the same relative phase). 

are in contact and the plate acceleration a(t) < —g. Fig- 
ure |^(a) compares the calculated flight time r ca i c with the 
measured flight time r eX pt as a function of T. (Figures 
H(c) and ||(d) indicate how r exp t is measured.) In the 
absence of parametric waves, the layer and ball motions 
are nearly the same: r oxp t = r ca i c . When parametric 
waves are present the layer is dilated at take-off, which 
reduces the effective take-off velocity of the layer and 
consequently decreases T exp t- Since 7 is calculated from 
model values of v c and tfn, it is only accurate when para- 
metric waves are absent (7 < Tent)- Because v c and t&t 
are both proportional to 1//, 7 is independent of /; this 
is consistent with Fig. || which shows that the parametric 
wave transitions depend only weakly on /. 

For r > 1, on each cycle the layer loses contact and 
later collides with the container. However, it is not until 
r w 2.4 that the parametric wave instability first occurs 
(7 = Tcrit) and the flat layer bifurcates to waves oscillat- 
ing at f/2, as shown in Fig. 1(a) p0| . Vertical lines ex- 
tending from Fig. ||(a) to Fig. 2 (b) mark the experimental 
values of T where transitions associated with paramet- 
ric waves occur. The f/2 wave patterns are squares for 
/ < 24 Hz and stripes for / > 40 Hz (Fig. g); at inter- 
mediate values of /, both patterns compete. For squares 
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the difference between the stability threshold for increas- 
ing and decreasing T is about twenty percent, while for 
stripes the hysteresis is small and perhaps even zero. 

Hexagonal patterns arise spontaneously from f/2 
waves when the vertical motion of the layer undergoes 
a period doubling bifurcation. Period doubling occurs 
when Texpt exceeds unity and subsequently becomes dou- 
ble valued (see Fig. ||(a)). As Fig. ||(e) shows, alternate 
collisions now occur before and after a c = —g; the layer 
is effectively forced at both / and f/2 because successive 
collisions occur with different v c and t^. Like f/2 waves, 
hexagonal patterns return to the same configuration after 
two periods but with an important difference: time trans- 
lation by one period is no longer equivalent to a spatial 
shift of one-half wavelength. Instead, two distinct pat- 
terns appear in alternate cycles — a set of isolated peaks 
on a triangular lattice become, on the next oscillation, 
hexagonal cells, each one centered on a former peak lo- 
cation. The amplitude of the cellular phase is maximum 
at the end of the long flight, while the amplitude of the 
peaked phase is maximum at the end of the short flight. 
Figure |2| shows that the T value for transition to hexagons 
(r w 3.9) is nearly constant for frequencies above where 
f/2 waves show a strong decrease in hysteresis. 

Spatial kinks separating domains of different relative 
phase appear at the period doubling bifurcation. Kinks 
arise because the period doubled motion is degenerate by 
7r in phase: a portion of the layer can begin its motion 
on one cycle or the next. For hexagons, this degeneracy 
manifests itself in the simultaneous appearance of both 
peaked and cellular spatial domains separated by a well- 
defined phase defect, as shown in Fig. |l|(b). 

As r is increased further, the amplitude of the hexag- 
onal pattern decreases until 7/7 cr it < 1 and parametric 
waves disappear. The layer now consists of flat domains 
connected by a kink, as shown in Fig. |l](c). Physically, 
this resonance corresponds to v c w 0: the layer lands 
gently on the plate. Flat domains with kinks arise solely 
from the period doubling instability and appear to be the 
2D analogue of the ID subharmonic instability in a glass 
bead layer reported by Douady et al. . 

Beyond onset of the flat domains with kinks, r cxpt is 
double valued as for the hexagonal patterns, but as T 
is increased the time of the shorter flight goes to zero 
and r exp t again becomes single valued. In both cases the 
motion is period doubled: in the former, flights are com- 
posed of long and short jumps (effective two-frequency 
forcing), while in the latter, flights can begin on one cy- 
cle or the next but the time between successive collisions 
is always 2/ f. In this latter state, T oxpt « T ca i c , indicating 
that the layer is compact at take-off. 

Parametric square and striped patterns reappear when 
again 7/7 cr it > 1 at T w 5.7. Because the time between 
collisions is two periods, the waves oscillate at //4 and, 
as Fig. [l](d) shows, appear in separate domains whose 
motions differ in phase by it (one domain is taking-off 



while the other is in mid-flight). Otherwise f/4 waves 
are the same as f/2 waves: the transition is accompanied 
by a decrease in r exp t with respect to r ca i c and squares 
appear subcritically at lower /, while at higher / stripes 
occur essentially without hysteresis. 

Although the layer motion is period doubled for f/4 
squares and stripes, hexagons do not appear because 
T Gxpt is single valued; the layer experiences the same 
conditions on consecutive collisions. However, when 
Texpt > 2, the layer motion undergoes a second period 
doubling — effective two frequency forcing (now at f/2 
and f/4) is restored and hexagons reappear. Figure |l|(e) 
shows that four different phases simultaneously exist in 
the system; each phase present before the second pe- 
riod doubling supports two separate phase domains of 
hexagons. As Fig. || indicates, the transition to f/4 
hexagons is non-hysteretic. 

The layer becomes spatially and temporally disordered 
for r « 7.6. Disorder is introduced by circular regions 
approximately three cells in diameter, which randomly 
appear in the hexagonal patterns and then shrink and 
disappear over the course of approximately 50//. The 
circular regions are 7r/4 out of phase with the domain 
that contains them and result from a single flight of tat ~ 
1//, instead of tat ~ 2// as for the rest of the domain. 
For r < 7.8 the four phase domains associated with the 
hexagonal pattern remain intact. However, for larger T, 
hexagons and the boundaries between the stable domains 
disappear; the layer consists of numerous small domains 
of wave-like structures with short spatial and temporal 
correlation (see Fig. |l](f)). The disordered state persists 
up to r = 14, the largest forcing examined. 

The ordered patterns we observe in granular media 
have similarities with those in vertically oscillated flu- 
ids (Faraday experiment) but an important differ- 
ence is that period doubling, which is not observed in 
fluids, leads to domains with different phase in the gran- 
ular layer. With this exception, it appears that the fac- 
tors determining pattern selection in the two media are 
similar even though, because of free flights and collisions, 
the effective forcing of the granular layer is different from 
the smooth forcing of the fluid. Dilation of the granular 
layer appears to be analogous to the inverse of viscosity 
in fluids because in the granular system we observe f/2 
squares for large dilation and stripes for small dilation, 
while in fluids squares are observed for small viscosity 
(less than 0.7 cm 2 /s) and stripes for large viscosity p2[ . 
In the Faraday experiment hexagons appear when the 
container is externally driven simultaneously at both / 
and f/2, but only for certain ranges of the relative phase 
and amplitude of the two drive signals; if the phase dif- 
ference is too small, squares are observed [Q . Hexagons 
appear in the granular layer when the intrinsic relative 
phase associated with the effective two-frequency forcing 
is large. In thinner particle layers the transition to ef- 
fective two frequency forcing occurs for smaller T; conse- 
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quently the relative phase between long and short flights 
is smaller and we observe patterns of squares with the 
same peak and cell structure as for hexagons. 

The pattern formation phenomena we have described 
result from the interaction of parametric waves and pe- 
riod doubling. In addition to the results reported here, 
we have performed experiments varying the particle resti- 
tution coefficient (0.5-0.95), density (2.3-11.4 g cm~ 3 ) 
and size (0.05-3 mm); the layer thickness (2-40 particles) 
and aspect ratio (2-100); and the pressure (10 _1 -10 -6 
Torr). No qualitative changes in the patterns or their 
associated bifurcations were observed. Moreover, in con- 
trast with most other cooperative dynamic phenomena in 
granular systems, we have found that transition thresh- 
olds and pattern organization become better defined as 
the pressure is reduced and the ratio of container size 
to pattern wavelength is increased; this indicates that 
interstitial gas and sidewall forcing do not organize the 
grain motion. Future work examining basic transport 
properties such as viscosity and momentum diffusivity in 
the robust patterns we have found should have practical 
consequences. 

The work of F. M. was supported by Fondecyt Grant 
No. 1941115. The work of P. B. U. and H. L. S. was 
supported by the Department of Energy, Office of Basic 
Energy Sciences Grant No. DE-FG03-93ER14312. The 
Texas-Chile collaboration was supported by a coopera- 
tive program between Conicyt and the National Science 
Foundation, Division of International Programs, Ameri- 
cas Program Grant No. INT-9415709. 



Mehta and J. M. Luck, Phys. Rev. Lett. 65, 393 (1990). 

[10] The parametric wave instability arises from two counter- 
propagating fluxes In. At onset in air with glass beads 
r ~ 3.3, while in vacuum F « 2.8. The strong increase 
in critical T in air near / = 40 Hz is not observed 
in vacuum. See B. Thomas, M. O. Mason, Y. A. Liu, 
and A. M. Squires, Powder Technol. 57, 267 (1989) for a 
description of other layer states observed in air. 

[11] S. Douady, S. Fauve, and C. Laroche, Europhys. Lett. 8, 
621 (1989). 

[12] W. S. Edwards and S. Fauve, J. Fluid Mech. 278, 123 
(1994). 

[13] H. W. Miiller, Phys. Rev. Lett. 71, 3287 (1993). 



TABLE I. Patterns and their associated instabilities. 
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